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Bosonic and fermionic Hubbard models on the checkerboard lattice are studied numerically for infinite on-site 
repulsion. At particle density n=l/4 and strong nearest-neighbor repulsion, insulating Valence Bond Crystals 
(VBC) of resonating particle pairs are stabilized. Their melting into superfluid/metallic phases under increasing 
hopping is investigated at T = QK. More specifically, we identify a novel and unconventional commensurate 
VBC supersolid region, precursor to the melting of the bosonic crystal. Hardcore bosons (spins) are compared 
to fermions (electrons), as well as positive to negative (frustrating) hoppings. 



Understanding the conditions and mechanisms for charge 
conduction in strongly correlated systems, especially frus- 
trated ones, is a highly non-trivial problem: aside from the 
strength of interactions, the lattice geometry, but also statistics 
of charge carriers, can play a relevant role. A case of interest is 
that of lattices with a pyrochlore-like structure: among those, 
various unusual electronic phases are encountered, from (in- 
sulating) spin liquids 1 1 1 to superconductivity, e.g. in spinels 
involving id electrons as quarter-filled LiTi204 |2|. Also, 
other spinels like LiV204 have heavy-fermion behaviors ||3| 
and are close to a Metal-Insulator (MI) transition f4l, raising 
attention to electronic correlations in a pyrochlore-based ge- 
ometry. 

On the other hand, properties of some half-polarized anti- 
ferromagnetic spinels |5 1 can be described similarly as those 
of anisotropic (XXZ) magnets on a pyrochlore lattice. In the 
simplified case of a spin- 1/2 XXZ magnet, representing up- 
(down-)spin by the presence (absence) of a boson, a map- 
ping to hard-core bosons with NN repulsive interactions can 
be used; an external magnetic field h^, analog of a bosonic 
chemical potential, tunes the bosonic density. Ath^ ^ Q and 
only an Ising coupling on the bonds of the checkerboard lat- 
tice of Fig. [T| (the two-dimensional analog of the pyrochlore 
lattice), the classical ground state (of density n=l/2) is highly 
degenerate and follows the so-called ice-rule constraint: the 
lowest Ising energy corresponds to precisely two bosons on 
every (flattened) tetrahedron. By tuning to reach a bosonic 
density of n=l/4, again the classical ground-states obey an ice- 
rule constraint of precisely one particle per tetrahedron (blue 
tetrahedra of Fig.[T]l. Second-order processes in the exchange 
coupling lead to the dynamics of a loop model ||6l|71 (n=l/2) or 
of a quantum dimer model (QDM) |8 1 (n=l/4). Note that sim- 
ilar lattice gas picture can also be relevant to easy-axis spin-S 
systems |9| (S* > 3/2). 

A specificity of these local constraints, aside from the finite 
ground state degeneracy at the classical level, is that charge 
excitations may a priori be fractionalized ifTOl : adding a par- 
ticle of charge — e to a system respecting initially the ice-rule 
creates 2 defects (tetrahedra on which this rule is violated) 
of fractional charge — e/2 which can move to neighboring 
tetrahedra (Fig. [l]). An insulating phase with fractional ex- 



citations could thus be considered, close enough to the MI 
transition such that fractionalization is allowed by kinetic en- 
ergy gain ifTTl . 
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FIG. 1: (Color online), a) Topological defects of fractional charge 
+e/2 (2 particles on a red tetrahedron) and —e/2 (empty blue tetra- 
hedron) and particle-hole excitation with two pairs of NN defects, 
b) Artistic view of the mixed- VBC supersolid phase at commensu- 
rate density n=l/4 showing inhomogeneous charge density (dots), 
resonating plaquettes (blue) and superfluid density (light blue). 

Lastly, we note that supersolidity has been predicted under 
doping a charge-ordered bosonic state fVT\. Commensurate 
supersolids can be observed as an out-of-equilibrium state of 
cold atoms |fT3l , or arising from the melting of a conventional 
charge-ordered insulator [H]. Recent developments of frus- 
trated optical lattices of cold atoms 1 15] open new directions. 

To address the above issues, in this Letter we consider Hub- 
bard hamiltonians on the checkerboard lattice in the infinite 
on-site repulsion (U) limit that reads: 



H=-tY, Pg{4,.c,,. + h.c.)PG + VY. 
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where cj ^ is a fermion (electron) creation operator at site i, 
a its spin index and ni is the onsite density operator. The 
Gutzwiller projector Pq enforces the U =oo limit (no doubly 
occupied site) 1 16*1. Note that the hopping term t and the near- 
est neighbor (NN) Coulomb repulsion V have the same mag- 
nitude on all bonds. We also consider hardcore bosons (spins 
- a bosonic density n < 1/2 corresponds to a spin system at 



2 



polarization 1 — 2n) replacing Ci^a by bosonic bi operators and 
dropping the sum on a. 

Although the lattice is frustrated, the t > bosonic model 
does not suffer from a sign problem and can be addressed us- 
ing Quantum Monte Carlo (QMC) techniques. In a closely 
related system, in which bosons hop on the square lattice 
and interact via bonds on the checkerboard lattice, analyti- 
cal arguments and finite-T QMC simulations have shown that 
the system (extrapolated to T = OK) undergoes, at density 
n=l/4, a weakly first-order phase transition from a superfluid 
(SF) phase at small V/t to a lai-ge V/t MI phase |17| (which 
breaks lattice translation symmetry), despite the possibility of 
being an unusual non-Landau transition flSl . A similar type 
of phase transition was also found on the Kagome lattice 1 19 |. 
Here, since the bosons can hop on all bonds of the lattice, 
we performed Green Function QMC (GFQMC) simulations 
directly at T = OK and in the canonical ensemble (n=l/4) 
on clusters up to 14a/2 x 14\/2 (N=392 sites) to determine 
the emerging quantum phases. For frustrated t < and/or 
the case of fermions, studies are restricted to smaller clusters, 
typically N=32 sites for a cluster compatible with all lattice 
symmetries) which can be handled using Lanczos Exact Diag- 
onalizations (ED) L20J . 
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FIG. 2: (Color online). Phase diagram of the t>Q bosonic hard-core 
Hubbard model as a function of the parameter t/V and the chemi- 
cal potential fi. Extrapolation to the thermodynamic limit (TDL) of 
T=OK GFQMC results in the canonical ensemble (dots) are com- 
pared to ED results. The two lobes correspond to the n=l/4 (bottom) 
and n=l/2 (top) VBC insulators (see text). A region of supersolid 
phase is found at commensurate n=l/4 filling (thick segment). 

A typical phase diagram of eq.|[T] is shown in Fig. [2] where 
the chemical potential /i has been introduced via a simple Leg- 
endre transformation w.r.t. the density. This phase diagram 
(for t>0 bosons, to be discussed later) raises a number of im- 
portant issues i.e. (i) the properties of the insulating VBC 
in the lobes, (ii) the nature of the phase transition towards 
the itinerant phases, (iii) the role of the quantum statistics 
(fermions vs bosons) and (iv) the effect of a frustrating t<0 
hopping. The paper is organized as follows. First, we investi- 
gate the evolution with increasing hopping of the gapped insu- 
lating phase towards a compressible fluid. Next, we character- 



ize the VBC insulator from various relevant structure factors 
and discuss the presence of a novel exotic supersolid phase 
with coexisting VBC and SF orders. Third, optics and dipolar 
excitations are discussed by calculating the optical properties 
in the VBC insulators and in the itinerant metallic (fermions) 
or SF (bosons) phases. Finally, we consider the physics of 
fractional defects beyond the ice-rule constraint. Note that the 
results shown below refer to ?i=l/4, i.e. to a quaiter-filled 
bosonic (1/8-filled fermionic) system. 

From the insulator to a compressible fiuid - The insulat- 
ing (incompressible) character of a system manifests itself 
by a finite gap in the quasiparticle charge excitation spec- 
trum. The gap closes in the metallic (fermions) or SF (bosons) 
phase fTTl|. Physically, the process of creating a particle-hole 
excitation as in Fig. [TJb) by moving a particle from one edge 
of the sample to the other edge costs a finite energy Ap, the 
charge gap, simply calculated by Ac = -E+ + — 2Eo with 
the (total) ground state (GS) energies Eq at filling n=l/4, and 
E± with one extra (less) particle. The ED results in Fig.[3|a) 
show that, both for bosons or fermions, Ac decreases linearly 
with \t\ from 2V at i = suggesting a transition at finite t 
to an itinerant GS. Despite the small cluster size some trends 
can be noticed; (i) the relative robustness of the insulator w.rt. 
its competing compressible phase does not seem to differ sig- 
nificantly between bosons and fermions and (ii) a frustrating 
t < hopping is much less efficient to destroy the insulator 
(see also Ref. 20). For unfrustrated bosons, GFQMC calcula- 
tion can be performed on much larger size clusters (up to 392 
sites), at T = OK and for very small ratio of t/V. A careful 
size-scaling enables us to extract the TDL of Ac, displayed in 
Fig.Qa). The transition characterized by the closure of Ac is 
located at t/V ~ 0.082, not so far from its estimation for the 
square hopping boson model iflTl . Note that the width A/i of 
the n=l/4 and n=l/2 insulating lobes of Fig.|2]is given by the 
extrapolated charge gap. 
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FIG. 3: (Color online). ED results obtained on a 32-site cluster 
for both fermions (empty symbols) and bosons (filled symbols): a) 
charge (circles) and optical (squares) gaps, all in units of V, b) order 
parameter M+{K), c) Drude weight and d) average number of — e/2 
charge defects. For fermions, discontinuities due to a level crossing 
at t/V ~ 0.035 might be a finite size effect. 
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VBC insulators and supersolid (SS) - As discussed above, 
at large V/t, mappings to generalized QDM describing the 
insulating phase are extremely useful to establish the VBC 
nature of the insulator, in particular its feature of resonat- 
ing particle pairs at n=l/4 ETl 1221 . For fermions, spin de- 
grees of freedom play an essential role, since the effective 
kinetic processes act only on singlet electron-pairs on the 
void plaquettes Il20ll22l . The insulators at n=l/2 and n=l/4 
have gapped triplet excitations and the same broken transla- 
tion symmetries 1221 as their bosonic VBC analogs of Refs.|6| 
and |2ll However, when approaching the phase transition by 
increasing t, the microscopic model ([T]l becomes necessary 
since the charge defects of Fig. [Tfa) play a crucial role. We 
then compute the structure factor of the diagonal operator 
P± = didj ± d^di defined in a previous work ll2Tl . where 
didj counts the pair (0 or 1) of particles facing each other 
across every void plaquette (on the same plaquette at 7r/2 an- 
gle for dkdi). Note that in the fermionic case, a projector 
on the 2-particle Sz=0 subspace is implicit in the definition 
of P±. Following the same symmetry considerations as in 
ETl . the s-wave (d-wave) structure factor at point K={tt,tt) 
(r=(0, 0)) signals plaquette ordering (rotational symmetry- 
breaking). The related order parameters M+{K) (plaquette) 
and M_ (F) (columnar) ll23l are defined as 



M±{q) 



(2) 



Data of Fig. |3jb) reveal striking differences between t > 
and t < with a weaker suppression of the VBC plaque- 
tte order parameter with \t\ in the frustrated case, following 
the same trend as Ac vs \t\. In the bosonic t > model, 
a careful size-scaling analysis of QMC data again enables to 
extract the TDL of the T=OK order parameters as shown in 
Fig.|4|b). This evidences clearly that these two order parame- 
ters of the mixed columnar-plaquette phase found in the large- 
V limit 121 J melt almost simultaneously at the phase transition 
around t/V ~ 0.089. Strikingly, the vanishing of the VBC or- 
der parameter appears at a larger value than the insulator-SF 
transition. This shows a persistence of the VBC order within 
the SF phase, i.e. a small region of SS phase. Note that (i) 
in contrast to the bosonic triangular lattice 1 12 1, this phase ap- 
pears here at a commensurate density, (ii) its VBC character 
strongly differs from usual charge ordering, and (iii) the pres- 
ence of kinetic processes on the crossed bonds is essential, 
the VBC insulator-SF transition being weakly first-order oth- 
erwise [17 1 with no intermediate phase. However the exact 
equality between the bond strengths within the crossed pla- 
quettes is not crucial and the exotic SS phase should exist in an 
extended parameter region. We speculate that, for fermions, a 
similar (narrow) exotic VBC metallic phase might also appear 
Optics - We now analyze the optical conductivity: first, it 
provides important information on transport, from its a; — > 
limit, the Drude weight (its finiteness characterizing a metal 
or SF). Note that for bosons at T=OK, this Drude weight cor- 
responds to the stiffness of superfluidity 1241 . Secondly, neu- 
tral dipolar excitations should appear in the finite frequency 
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FIG. 4: (Color online). GFQMC results for (hard-core) bosons with 
t>0. TDL of the charge gap in units of V (a), the order parameters 
M+{K) and A'/_(r) with the Brillouin zone depicted in the inset 
(b), the stiffness parameter rescaled by the number of sites TV (c) and 
the average number of — e/2 charge defects (d). The same scale is 
used to plot these quantities versus t/V to emphasize the transition 
point close to 0.085. Size-scalings are depicted in the insets. 



absorption spectrum. Observe that, for fermions, these excita- 
tions are also spin singlets. First, we calculate by Lanczos ED 
the finite frequency conductivity. 



uj ^ uj + iO+ — H + Lq ' 



1 



where is the current operator along one of the diagonal di- 
rections as shown on Fig.[T] The data on Fig. [5] clearly show, 
at small an optical gap Aopt ^ V corresponding to an 

absorption energy threshold. As shown in Fig. [3ja), this gap 
follows closely - for small enough \t\/V - the linear behav- 
ior with |i| of Ac. We beUeve the up- turn at t/y 0.08 
indicates the transition to the metallic/superconducting phase, 
where the energy scale of excitations cross-over from V to t. 

Interestingly, low-energy peaks are observed in a{uj) below 
Aopt 31 energies /V set by the effective QDM hamilto- 
nian. While such q={0, 0) low-energy dipolar excitations cor- 
respond to usual excitations of a metal/SF, it is not clear yet 
whether they would also survive in the TDL in the insulator, 
as occuring from local charge fluctuations. 

The Drude weight Ds, extracted from the above ED re- 
sults via the optical sum rule and plotted in Fig.|3jc), increases 
rapidly with \t\/V but does not allow to locate the insulator- 
metal (or SF) transition for such small systems. Fortunately, 
in the case of bosons and t>0, it is possible to compute the 
Drude weight with the GFQMC method using the winding 
numbers of the flow of the particles in imaginary time t |24|. 
At a continuous solid to SF transition, the stiffness scales as 
Ds = L-^'fiL^/" {t/V - {t/V)c) , tL^) with the dynami- 
cal z and the correlation length v exponents 1 19 1. Searching 
for a universal behavior of / w.rt. L, we obtain a very accu- 
rate estimation of the transition at t/V ~ 0.082 (crossing in 
Fig.|4|c)), with z = 2 and = 1/2. In agreement with the 
vanishing of Ac, this is an extra evidence in favor of a win- 
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FIG. 5: (Color online). (j{uj) for a) the fermions and b) the bosons 
at different values of t/V (ED on a 32-site cluster). Insets: relative 
weight of the low-energy dipolar matrix elements (only visible for 
fermions on the a(uj) plot), plotted vs t/V in log-log scale. 

dow of width A(i/V^) ~ 7 • 10"^ of a SS phase exhibiting a 
finite SF stiffness Ds and = but with two finite VBC 
order parameters. 

The physics of fractional defrcts - Within the effective 
QDM, the ampHtude of dimer flips {i.e. simuhaneous 2- 
particle move on a void plaquette |6|) goes like /V . Hence, 
the concentration of ±e/2 defects (c/Fig. [T|a)) shows in 
Fig. [3] for very smafl t, a {t/VY behavior, symmetrically for 
t > and i < 0. However, as soon as \t\/V > 0.03, the 
observed behaviors are beyond the physics of the effective 
QDM. Nevertheless, the phase transition occurs at a small t/y 
value so that the ice-rule constraint should still operate in the 
metallic/SF phase which contains only ^6-10% of fractional 
defects. Since the deconfinement |17| of these defects occurs 
at the vanishing of the VBC order parameter, they should still 
be confined in the SS phase. 

We summarize here the main findings of this work. First, 
a strong asymmetry is found between t> and t<0, a frus- 
trated hopping (t<0) being less efficient to melt the VBC 
insulator. This originates from a relatively smaller increase 
with \t\ of the average kinetic energy and concentration no 
of fractional charged defects, as shown in Fig. |3jd). Sec- 
ondly, for the non frustrated t>0 bosonic model, evidence 
are provided for a double quantum phase transition between 
the VBC insulator at small t and the SF at large t, involving 
an intermediate commensurate f25l supersolid phase breaking 
both translation and rotational symmetries. Its VBC charac- 
ter (anisotropic resonating plaquettes) clearly distinguishes it 
from its charge ordered counterpart |14|. We speculate that 
such a feature should be present in the frustrated t<0 case 
and/or for fermions (with an intermediate VBC metal). Lastly, 
optical properties are investigated, showing an optical gap fol- 
lowing the behavior of the charge gap. 
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